Structural identifiability of dynamic systems biology models

22 páginas, 5 figuras, 2 tablas.-- This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.A powerful way of gaining insight into biological systems is by creating a nonlinear differential equation model, which usually contains many unknown parameters. Such a model is called structurally identifiable if it is possible to determine the values of its parameters from measurements of the model outputs. Structural identifiability is a prerequisite for parameter estimation, and should be assessed before exploiting a model. However, this analysis is seldom performed due to the high computational cost involved in the necessary symbolic calculations, which quickly becomes prohibitive as the problem size increases. In this paper we show how to analyse the structural identifiability of a very general class of nonlinear models by extending methods originally developed for studying observability. We present results about models whose identifiability had not been previously determined, report unidentifiabilities that had not been found before, and show how to modify those unidentifiable models to make them identifiable. This method helps prevent problems caused by lack of identifiability analysis, which can compromise the success of tasks such as experiment design, parameter estimation, and model-based optimization. The procedure is called STRIKE-GOLDD (STRuctural Identifiability taKen as Extended-Generalized Observability with Lie Derivatives and Decomposition), and it is implemented in a MATLAB toolbox which is available as open source software. The broad applicability of this approach facilitates the analysis of the increasingly complex models used in systems biology and other areasAFV acknowledges funding from the Galician government (Xunta de Galiza, Consellería de Cultura, Educación e Ordenación Universitaria http://www.edu.xunta.es/portal/taxonomy/term/206) through the I2C postdoctoral program, fellowship ED481B2014/133-0. AB and AFV were partially supported by grant DPI2013-47100-C2-2-P from the Spanish Ministry of Economy and Competitiveness (MINECO). AFV acknowledges additional funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No 686282 (CanPathPro). AP was partially supported through EPSRC projects EP/M002454/1 and EP/J012041/1.Peer reviewe

Data-driven reverse engineering of signaling pathways using ensembles of dynamic models

Signaling pathways play a key role in complex diseases such as cancer, for which the development of novel therapies is a difficult, expensive and laborious task. Computational models that can predict the effect of a new combination of drugs without having to test it experimentally can help in accelerating this process. In particular, network-based dynamic models of these pathways hold promise to both understand and predict the effect of therapeutics. However, their use is currently hampered by limitations in our knowledge of the underlying biochemistry, as well as in the experimental and computational technologies used for calibrating the models. Thus, the results from such models need to be carefully interpreted and used in order to avoid biased predictions. Here we present a procedure that deals with this uncertainty by using experimental data to build an ensemble of dynamic models. The method incorporates steps to reduce overfitting and maximize predictive capability. We find that by combining the outputs of individual models in an ensemble it is possible to obtain a more robust prediction. We report results obtained with this method, which we call SELDOM (enSEmbLe of Dynamic lOgic-based Models), showing that it improves the predictions previously reported for several challenging problems.JRB and DH acknowledge funding from the EU FP7 project NICHE (ITN Grant number 289384). JRB acknowledges funding from the Spanish MINECO project SYNBIOFACTORY (grant number DPI2014-55276-C5-2-R). AFV acknowledges funding from the Galician government (Xunta de Galiza) through the I2C postdoctoral fellowship ED481B2014/133-0. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.info:eu-repo/semantics/publishedVersio

An Integrated Qualitative and Quantitative Biochemical Model Learning Framework Using Evolutionary Strategy and Simulated Annealing

The authors would like to thank the support on this research by the CRISP Project (Combinatorial Responses In Stress Pathways) funded by the BBSRC (BB/F00513X/1) under the Systems Approaches to Biological Research (SABR) Initiative.Peer reviewedPublisher PD

Fire and brief human occupations in Iberia during MIS 4: Evidence from Abric del Pastor (Alcoy, Spain)

There is a relatively low amount of Middle Paleolithic sites in Europe dating to MIS 4. Of the few that exist, several of them lack evidence for anthropogenic fire, raising the question of how this period of global cooling may have affected the Neanderthal population. The Iberian Peninsula is a key area to explore this issue, as it has been considered as a glacial refugium during critical periods of the Neanderthal timeline and might therefore yield archaeological contexts in which we can explore possible changes in the behaviour and settlement patterns of Neanderthal groups during MIS 4. Here we report recent data from Abric del Pastor, a small rock shelter in Alcoy (Alicante, Spain) with a stratified deposit containing Middle Palaeolithic remains. We present absolute dates that frame the sequence within MIS 4 and multi-proxy geoarchaeological evidence of in situ anthropogenic fire, including microscopic evidence of in situ combustion residues and thermally altered sediment. We also present archaeostratigraphic evidence of recurrent, functionally diverse, brief human occupation of the rock shelter. Our results suggest that Neanderthals occupied the Central Mediterranean coast of the Iberian Peninsula during MIS 4, that these Neanderthals were not undergoing climatic stress and they were habitual fire users.This research was funded by a Leakey Foundation General Grant, Spanish Ministry of Science, Innovation and Universities Projects HAR2008-06117/HIST and HAR2015-68321-P, Junta de Castilla y León-FEDER Project BU235P18, the LabEx Sciences Archéologiques de Bordeaux (LaScArBx ANR-10-LABX-52) and ERC Consolidator Grant ERC-CoG-2014. Archaeological excavations at Abric del Pastor are supported by the Archaeological Museum of Alcoy and the Government of Valencia Cultural Heritage Department

Non-CO2 generating energy shares in the world : cross-country differences and polarization.

The aim of this paper is to examine the spatial distribution of non-CO2 generating energy sources in the world for the period 1990–2009, paying special attention to the evolution of cross-country disparities. To this end, after carrying out a classical convergence analysis, a more thorough investigation of the entire distribution is presented by examining its external shape, the intra-distribution dynamics and the long-run equilibrium distribution. This analysis reveals the existence of a weak, rather insignificant, convergence process and that large crosscountry differences are likely to persist in the long-run. Next, as polarization indicators are a proper way of appraising potential conflict in international environmental negotiations, we test whether, or not, the distribution dynamics concurs with the presence of polarization. Our results indicate that two poles can be clearly differentiated, one with high and other with low non-CO2 generating energy shares. In view of these findings, and to ensure a fair transition to a sustainable energy system, the paper calls for the development of an ambitious clean energy agenda, especially in countries with low non-CO2 generating energy shares

Multiplexing information flow through dynamic signalling systems

We consider how a signalling system can act as an information hub by multiplexing information arising from multiple signals. We formally define multiplexing, mathematically characterise which systems can multiplex and how well they can do it. While the results of this paper are theoretical, to motivate the idea of multiplexing, we provide experimental evidence that tentatively suggests that the NF-κB transcription factor can multiplex information about changes in multiple signals. We believe that our theoretical results may resolve the apparent paradox of how a system like NF-κB that regulates cell fate and inflammatory signalling in response to diverse stimuli can appear to have the low information carrying capacity suggested by recent studies on scalar signals. In carrying out our study, we introduce new methods for the analysis of large, nonlinear stochastic dynamic models, and develop computational algorithms that facilitate the calculation of fundamental constructs of information theory such as Kullback–Leibler divergences and sensitivity matrices, and link these methods to a new theory about multiplexing information. We show that many current models such as those of the NF-κB system cannot multiplex effectively and provide models that overcome this limitation using post-transcriptional modifications

Early holocenic and historic mtDNA african signatures in the iberian peninsula: The andalusian region as a paradigm

Determining the timing, identity and direction of migrations in the Mediterranean Basin, the role of "migratory routes" in and among regions of Africa, Europe and Asia, and the effects of sex-specific behaviors of population movements have important implications for our understanding of the present human genetic diversity. A crucial component of the Mediterranean world is its westernmost region. Clear features of transcontinental ancient contacts between North African and Iberian populations surrounding the maritime region of Gibraltar Strait have been identified from archeological data. The attempt to discern origin and dates of migration between close geographically related regions has been a challenge in the field of uniparental-based population genetics. Mitochondrial DNA (mtDNA) studies have been focused on surveying the H1, H3 and V lineages when trying to ascertain north-south migrations, and U6 and L in the opposite direction, assuming that those lineages are good proxies for the ancestry of each side of the Mediterranean. To this end, in the present work we have screened entire mtDNA sequences belonging to U6, M1 and L haplogroups in Andalusians--from Huelva and Granada provinces--and Moroccan Berbers. We present here pioneer data and interpretations on the role of NW Africa and the Iberian Peninsula regarding the time of origin, number of founders and expansion directions of these specific markers. The estimated entrance of the North African U6 lineages into Iberia at 10 ky correlates well with other L African clades, indicating that U6 and some L lineages moved together from Africa to Iberia in the Early Holocene. Still, founder analysis highlights that the high sharing of lineages between North Africa and Iberia results from a complex process continued through time, impairing simplistic interpretations. In particular, our work supports the existence of an ancient, frequently denied, bridge connecting the Maghreb and Andalusia.Financial support was provided by the Spanish Ministry of Competitiveness through Research Project CGL2010-15191/BOS granted to RC and International Mobility Program Acciones Integradas Hispano-Portuguesas (PRI-AIBPT-2011-1004) granted to RC (Spain) and LP (Portugal) (http://www.mineco.gob.es/portal/site/mineco/idi). The E.C. Sixth Framework Programme under Contract n° ERAS-CT-2003-980409 (EUROCORES project of the European Science Foundation) also provided financial support to JMD for North African population research. CLH has a predoctoral fellowship granted by Complutense University. PS is supported by FCT Investigator Programme (IF/01641/2013). IPATIMUP (https://www.ipatimup.pt/) integrates the Instituto the Investigação em Saúde (i3S) Research Unit, which is partially supported by FCT, the Portuguese Foundation for Science and Technology. IPATIMUP is funded by FEDER funds through the Operational Programme for Competitiveness Factors - COMPETE and National Funds through the FCT - under the project PEst-C/SAU/LA0003/2013. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript

Cattle genome-wide analysis reveals genetic signatures in trypanotolerant N'Dama

Abstract Background Indigenous cattle in Africa have adapted to various local environments to acquire superior phenotypes that enhance their survival under harsh conditions. While many studies investigated the adaptation of overall African cattle, genetic characteristics of each breed have been poorly studied. Results We performed the comparative genome-wide analysis to assess evidence for subspeciation within species at the genetic level in trypanotolerant N’Dama cattle. We analysed genetic variation patterns in N’Dama from the genomes of 101 cattle breeds including 48 samples of five indigenous African cattle breeds and 53 samples of various commercial breeds. Analysis of SNP variances between cattle breeds using wMI, XP-CLR, and XP-EHH detected genes containing N’Dama-specific genetic variants and their potential associations. Functional annotation analysis revealed that these genes are associated with ossification, neurological and immune system. Particularly, the genes involved in bone formation indicate that local adaptation of N’Dama may engage in skeletal growth as well as immune systems. Conclusions Our results imply that N’Dama might have acquired distinct genotypes associated with growth and regulation of regional diseases including trypanosomiasis. Moreover, this study offers significant insights into identifying genetic signatures for natural and artificial selection of diverse African cattle breeds

Structural Identifiability Analysis via Extended Observability and Decomposition

Structural identifiability analysis of nonlinear dynamic models requires symbolic manipulations, whose computational cost rises very fast with problem size. This hampers the application of these techniques to the large models which are increasingly common in systems biology. Here we present a method to assess parametric identifiability based on the framework of nonlinear observability. Essentially, our method considers model parameters as particular cases of state variables with zero dynamics, and evaluates structural identifiability by calculating the rank of a generalized observability-identifiability matrix. If a model is unidentifiable as a whole, the method determines the identifiability of its individual parameters. For models whose size or complexity prevents the direct application of this procedure, an optimization approach is used to decompose them into tractable subsystems. We demonstrate the feasibility of this approach by applying it to three well-known case studies

Structural Identifiability Analysis via Extended Observability and Decomposition

Structural identifiability analysis of nonlinear dynamic models requires symbolic manipulations, whose computational cost rises very fast with problem size. This hampers the application of these techniques to the large models which are increasingly common in systems biology. Here we present a method to assess parametric identifiability based on the framework of nonlinear observability. Essentially, our method considers model parameters as particular cases of state variables with zero dynamics, and evaluates structural identifiability by calculating the rank of a generalized observability-identifiability matrix. If a model is unidentifiable as a whole, the method determines the identifiability of its individual parameters. For models whose size or complexity prevents the direct application of this procedure, an optimization approach is used to decompose them into tractable subsystems. We demonstrate the feasibility of this approach by applying it to three well-known case studies